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Abstract 



A search for the decay Ks — > tt^ fx^ ix~ has been made by the NA48/1 Collabora- 
tion at the CERN SPS accelerator. The data were collected during 2002 with a 
high-intensity Ks beam. Six events were found with a background expectation of 
0.22^q}^ events. Using a vector matrix element and unit form factor, the measured 
branching ratio is 

B{Ks ttV+M") = [2.9l|J(stoi) ± Q.2{syst)\ x 10"^. 



1 Introduction 

This paper reports the first observation of the decay Kg — > 7r°//'^/i~ and a measure- 
ment of its branching ratio. The analysis was carried out on data taken during 2002 by the 
NA48/1 experiment at the CERN SPS, which also recently reported the first observation 
of the decay Ks — > 7r°e"'"e~ [1] . 

The physics interest of the Ks — > 7r'^/x+/x^ decay is that it measures the indirect 
CP violating contribution of the decay Kl — > /i'^ , thereby allowing the direct CP 
violating component of the Kl decay to be extracted. This can provide input to the de- 
termination of the imaginary part of the element Vtd of the Cabibbo, Kobayashi, Maskawa 
(CKM) mixing matrix. In addition, the decay Ks —>■ Ti^ii^ ii~ can be used to study the 
structure of the K — > 777* form factor. This is particularly interesting if combined with 
the Ks — > 7r°e"'"e~ results since both decays are expected to be dominated by the exchange 
of a single virtual photon [K — > 777* — > 7r£+£~). 

2 Experimental Setup 

The beam line and detector built by the NA48 Collaboration to measure the Re{e' /e) 
parameter [2] from a comparison of Ks,l tt+tt^, tt^tt" decays was used, with the modi- 
fications described below. 

2.1 Beam 

The experiment was performed at the CERN SPS accelerator and used a 400 GeV 
proton beam impinging on a Be target to produce a neutral beam. The spill length was 
4.8 s out of a 16.8 s cycle time. The proton intensity was fairly constant during the spill 
with a mean of 5 x 10^° protons per spill. 

The neutral kaon beam line was modified as follows: the Kl beam was blocked and 
a small, additional sweeping magnet was installed above the Ks collimator. In order to 
reduce the number of photons in the neutral beam, primarily from 7r° decays, a platinum 
absorber 24 mm thick was placed in the beam between the target and the main sweeping 
magnet which deflected charged particles into the collimator. The 5.1m long Ks collima- 
tor, with its axis at an angle of 4.2 mrad to the proton beam direction, selected a beam of 
neutral long-lived particles {Ks, Kl, A", n, 7's, etc.). On average, there were 2 x 10^ 
Ks decays per spill in the fiducial volume downstream of the collimator, with momenta 
between 60 and 200 GeV (and a mean energy of ~ 110 GeV). An average of 0.017 Kl 
decays were expected for every Ks decay within the first three Ks lifetimes from the end 
of the collimator. 



2.2 Detector 

In order to minimize interactions of the neutral beam, the coUimator was imme- 
diately followed by a 90 m long evacuated tank. The tank was terminated by a 0.3% 

radiation length (Xq) thick Kevlar window, except in a region close to the beam which 
continued in a vacuum pipe through the centre of the downstream detectors. 

2.2.1 Tracking 

The tracking was performed with a spectrometer housed in a helium gas volume. 
It consisted of two drift chambers before and two after a dipole magnet with a horizontal 
transverse momentum kick of 265 MeV. Each chamber had four views, each of which 
had two sense wire planes. The resulting space points were typically reconstructed with 
a resolution of 150 yum in each projection. The spectrometer momentum resolution was 
parameterised as: 

ap/p = 0.48% ® 0.015% • p, 

where p is in GeV. This gave a resolution of 3 MeV on the reconstructed kaon mass in 
Ks — > 7r+7r~ decays. The track time resolution was 1.4 ns. 

2.2.2 Electromagnetic Calorimetry 

The detection and measurement of electromagnetic showers were performed with a 
27 Xq deep liquid krypton calorimeter (LKr). The energy resolution was parameterised as 
[3]: 

a(E)/E = 3.2%/ ® 9%/E ® 0.42%, 

where E is in GeV. The calorimeter was subdivided into 13,500 cells of transverse dimen- 
sion 2 cm X 2 cm, which resulted in a transverse position resolution better than 1.3 mm 
for single photons with energy above 20 GeV. The 7r° mass resolution was 0.8 MeV, while 
the time resolution of the calorimeter for a single shower was better than 300 ps. 

2.2.3 Scintillator Detectors and Muon Detector 

A scintillator hodoscope was located between the spectrometer and the LKr. It 
consisted of two planes, segmented in horizontal and vertical strips respectively, with 
each plane arranged in four quadrants. The time resolution for the hodoscope system 
was 200 ps. Downstream of the LKr calorimeter was an iron-scintillator sandwich hadron 
calorimeter (HCAL), followed by muon counters (MUG) which consisted of three planes of 
plastic scintillators, each shielded by an 80 cm thick iron wall. The first two planes, MIX 
and MIY, consisted of 25 cm wide horizontal [MIX) and vertical [MIY) scintillator 
strips, with a length of 2.7 m. The third plane, M2X , consisted of horizontal strips of 
width 44.6 cm and was mainly used to measure the efficiency of the MIX and MIY 
counters. The central strip in each plane was split with a gap of 21 cm to accommodate 
the beam pipe. The fiducial volume of the experiment was principally determined by 
the LKr calorimeter acceptance, together with seven rings of scintillation counters which 
surrounded the decay volume to veto activity outside this region. 

2.2.4 Trigger and Readout 

The detector was sampled every 25 ns, and samples were recorded in time windows of 
200 ns. The extended time window allowed the rate of accidental activity to be investigated 
from the time sidebands. 



The trigger selection for Ks v: ji'^ candidates consisted of a first-stage hard- 
ware trigger followed by a second-stage software trigger. The hardware trigger, for the 
first 40% of the data taken, selected events satisfying the following conditions: 

— at least one hit in each of the horizontal and vertical planes of the hodoscope, within 

the same quadrant; 

— at least three hit wires in at least three views of DCHl integrated over 200 ns; 

— a track vertex located within 90 m of the end of the collimator; 

— no hit in the two ring scintillator counters farthest downstream; 

— a two-muon signal (2/itight) from the muon counters. 

The 2/itight signal required at least two hits in each of the first two muon counter planes 
{MIX and MIY). 

For the last 60% of the data taken, a second, parallel, hardware trigger component 
was added in which the 2//tight condition was replaced by the following requirements: 

— hadron calorimeter energy less than lOGeV; 

— electromagnetic calorimeter energy greater than 15 GeV; 

— a two-muon signal (2/iioosc) from the muon counters. 

The 2/iioose signal was similar to the 2/itight signal but allowed one of the first two muon 
planes [MIX or MIY) to contain only a single hit. The addition of the 2/tioose trigger 
significantly improved the trigger acceptance for the Kg 7r°/t+/(~ signal, as discussed 
below in Section 3.4. 

The software trigger (for the entire data set) required: 

— at least two tracks in the drift chambers that were not associated with ener- 
getic (>5 GeV) clusters in the LKr; 

— less than 10 GeV in the hadronic calorimeter; 

— at least one hit in the first two muon counter planes. 

Events that satisfied the trigger conditions were recorded and reprocessed with improved 
calibrations to obtain the final data sample. 

3 Data Analysis 

The signal channel Ks 7r°/i"'"/i~ required the identification of two muons of op- 
posite charge accompanied by two additional clusters (photons) in the LKr. The kaon 
energy was estimated from the sum of the track momenta and cluster energies and was 
required to be between 60 and 200 GeV. 

The invariant mass, m^^, for the two muons from a Kg — > 7r°/( decay is limited 
by kinematics to: 

0.211 GeV < m^^ < 0.363 GeV 

where the lower limit corresponds to m^-\-m^, while the upper limit is given by mft- — m^o. 

The known Kg mass was used as a constraint (sec Eq. 1 below) to reconstruct the 
longitudinal vertex position using information from the charged tracks and the clusters. 
The invariant mass, m-y-y, of the two photons was reconstructed using the measured posi- 
tions and energies of the two clusters, assuming this vertex position. The total invariant 
mass, m^^TT, was calculated using the vertex reconstructed from the charged tracks alone, 
with the known 7r° mass imposed as a constraint to improve the kaon mass resolution. 

3.1 Signal and Control regions 

The similarity between the kinematics of the Kl — > tt^tt+tt" and Ks tt^ii'^ij,' 
channels allowed the former channel to be used to evaluate the resolutions (7m„_ and (7m„,_ 



on m^-y and m^^^, respectively. These were measured to be cr^^^ = 0.8 MeV and cr^^^^ = 
3 MeV, in agreement with a detailed Monte Carlo simulation based on GEANT [4] . 
A signal region and a control region were defined in the (m-y-y, m^^^) plane: 

- SIGNAL REGION: 

\m^^ - m^o| < 2.5 (T™^^ and |m^^^ - m^l < 2.5 cr^^^^; 

- CONTROL REGION: 

3cr„^^^ < \m^^ - m^o| < Qa^^^ and 3cr^^^^ < |m^^^ - mK\ < 6cr^^^„. 
The signal and control regions were kept masked until the analysis cuts were finalised in 
order to keep the event selection unbiased. The cuts to reject the background were set 
using both data and Monte Carlo simulation. An extended control region was used to 
estimate the contributions from various backgrounds. 



3.2 Event selection 

In order to identify muons, charged tracks reconstructed in the spectrometer were 
extrapolated to the MUC planes and associated with MUC hits using cuts on the spatial 
separation and time difference between the extrapolated track and the MUC hit. Multiple 
scattering was taken into account before applying the spatial cut, and light propagation 
along the MUC strips was taken into account before applying the time-difference cut. 
For muons with momentum greater than 10 GeV, the efficiency of the MUC counters was 
found to be 0.99 ± 0.01. 

Track pairs reconstructed in the spectrometer were classified as muons if the fol- 
lowing conditions were satisfied: 

— either no LKr cluster was associated to each track, or the energy of any associated 
cluster was less than 1.5 GeV. A cluster was associated with a track if the separation 
of the cluster and the track was less than 6 cm and the absolute time difference 
between the track time and the cluster time was less than 4.5 ns; 

— each track had associated hits in the first two planes {MIX and MIY) of the MUC; 

— the total energy in the hadron calorimeter was less than 10 GeV; 

— for each track, the time measured by either the drift chambers or the trigger ho- 
doscope and the time measured by the muon detector did not differ by more than 
4.5 ns. In about 95% of the events the track time was given by the trigger hodoscope, 
while for the remaining events the time was given by the drift chambers; 

— the average of the track times and the average time in the muon detector did not 
differ by more than 3 ns; 

— the tracks had opposite charge; 

— the individual track momenta were above 10 GeV in order to have high muon 
identification efficiency; 

— the transverse momentum with respect to the kaon line of fiight of each track was 
below 0.180 GeV; 

— the distance between the impact points of the two tracks at the first drift chamber 
was greater than 10 cm in order to have high trigger vertex reconstruction efficiency; 

— the ratio of the track momenta was larger than 1/3 or smaller than 3 in order to 
reduce the background from A and A decays coming either directly from the target 
or from — * A(A)7r° decays; 

— the reconstructed invariant mass assuming A —>■ p7T~ (or A) was not compatible 
with the A(A) mass ( 1.1115 GeV < rup^ < 1.1200 GeV); 

— 0.21 GeV < m„„ < 0.36 GeV. 



Clusters reconstructed in the LKr calorimeter were classified as photons if the following 
conditions were satisfied: 

— the energy of the cluster was greater than 3 GeV and less than 100 GeV; 

— the distance to the nearest cluster was greater than 10 cm in order to minimize the 
effect of energy sharing on cluster reconstruction; 

— the distance to the nearest extrapolated track was greater than 20 cm. 
In addition, a 7r° candidate satisfied the following conditions: 

— the total energy of the clusters was greater than 20 GeV; 

— the transverse momentum with respect to the kaon line of fiight of the 7r° candidate 
was between 0.05 and 0.25 GeV. 

Three quantities related to the decay vertex were computed: 

— total vertex, zk- 

To compute the longitudinal vertex position, zk, the kaon mass was assumed and 
the kinematical information from all the particles involved in the decay was used. 
The distance from the vertex to the LKr calorimeter, (Ik-i was calculated as: 



, ^ / T.^^jPiPid% 

^ \l ml -2ml -2 T.^^o^E^E^ ' PiPj) 
where dij is the separation between the i'th and j'th particles in the LKr plane after 
linearly extrapolating charged tracks from before the spectrometer magnet, and Ei 
and Pi are the energy and the momentum of the i'th particle. The vertex is given 
i>J Zk = z (LKr position) — dx- This formula follows from energy and momentum 
conservation applied to the decay K'^ — > /I'^ii'^^. The approximation cos a 



ij 



1 — I {dij/dK^ was made, where aij is the angle defined between the momenta of 
each pair of particles. 

— 7r° vertex, z-^^o . 

The 7r° vertex position along the beam direction, z^^o, was computed in a similar 
way to the total vertex, but used only the two photon clusters and imposed the 7r° 
mass instead of the kaon mass as a constraint. 

— track vertex, Zch- 

The track vertex was computed by finding the position of the closest distance of 
approach (CDA) between the two tracks. The CDA was required to be less than 
1.5 cm in order to reject muons from pion decay. 
Further cuts were applied: 

— on the assumption that the observed event was a kaon decay, its proper lifetime was 
calculated using the zk, Zt^o and Zch vertices, taking the end of the final collimator 
as the origin. All three calculated proper lifetimes were required to be between 
and 3 Kg mean lifetimes. 

— the energy- weighted centre-of-gravity (COG) of the event was required to be less 
than 5 cm in order to further reduce the background from accidental activity. The 
COG was defined as COG=-^ (^^ E^xiY + (Si ^iVi)'^/ Si where Xi and yi are 
the coordinates of the i'th particle in the LKr plane after linearly extrapolating 
charged tracks from before the spectrometer magnet, and is the energy of the 
i'th particle. 

— events with extra tracks within the trigger window or extra clusters in the readout 
window were rejected. 

— the time difference At between the average track time and the average time of the 
LKr clusters forming the tt*^ was required to be less than 1.5 ns. 



— fiducial cuts: tracks were required to be at least 12 cm from the centre of the drift 
chambers and clusters in the LKr were required to be at least 15 cm from the centre 
of the LKr, 11cm from its outside borders, and 2 cm from any defective cell. 

After applying all these cuts, the signal and control regions were unmasked. Six events 
were found in the signal region, and none in the control region (see Fig. 1). The kinematical 
quantities for these events are summarised in Table 1. 

3.3 Background 

There are two categories of background: physical background and accidental back- 
ground. Physical background was defined as that due to a single kaon decay, for example 
Kl — >■ 7r"'"7r~7r°, where the two charged pions decayed into muons. Accidental background 
was caused by the overlap of particles coming from two separate decays that happened to 
be in time and faked the signal, for example, an overlap between the decays Kl ttiii/, 
where vr — > fiu, and Ks tt^tt^, where two photons missed the detector. The accidental 
background contained two components: (a) overlapping fragments from different proton 
interactions in the target; and (b) fragments from associated production [KK or AK) due 
to a single proton interaction. The first component could be estimated through a study 
of out-of-time sidebands. 

3.3.1 Physical backgrounds 

Only two sources of physical background were found to give a significant contribu- 
tion: 

(1) Kl,s l^'^l^'ll decays: 

The Ki — > /i"'"/i~77 branching ratio has been measured by the KTeV Collaboration to 
be (l.O^o'g) X 10~^ [5]. Using Monte Carlo simulation, the acceptance for Kl — > /i"'";U~77 
decays in NA48 was found to be 5 x 10~^. This small value originates from a geometrical 
acceptance of about 20% together with the low probability that the invariant mass of the 
two photons is consistent with the 7r° mass. The Kg decay into the same final state was 
taken into account by assuming equal decay rates for Ks and Kl. The total background 
was estimated to be 0.04+[!:j^3 event. 

(2) Kl 7r"'"7r~7r° decay, where both pions have decayed in fiight: 

For this background, the reconstructed total mass m^^^ lies below the known kaon mass 
as a consequence of the missing energy of the neutrinos and the use of the muon mass 
rather than the pion mass in the reconstruction. For the same reasons, the z vertex 
position (Eq. 1) computed by imposing the K^ mass constraint is pushed upstream from 
the actual vertex position of the tt+tt'tt^ decay. The shift in the z vertex position leads 
to a reconstructed two-photon mass m^j lying above the known 7r° mass. These effects 
are clearly seen in Fig. 2, where the m^^ and m^^^r distributions are compared for data 
and Monte Carlo, after applying all cuts and after removing the Ks proper lifetime cuts 
calculated from zk- The Monte Carlo gives a good description of the data. 

An extended control region was defined in the {ninnT^jm^^) plane, between 6a and 
12(7 around the 7r° and Ks masses. No data event was found in this region. A Monte 
Carlo sample equivalent to 24 times the 2002 data was used to estimate the background. 
The cuts on the Ks lifetime, on the transverse momenta of the muons and the pion and 
on the total energy of the two photons were relaxed, giving 11 Monte Carlo events in 
the extended control region and no event in the signal region. This led to a background 
estimate of < 2.44/24 = 0.10 at 90%CL. It was checked that the transverse momentum 
and photon energy cuts did not significantly affect the shape of the relevant background 



distributions. After re-applying all cuts except the Ks lifetime cut, 2 of the 11 Monte Carlo 
events remained in the extended control region. This led to an additional extrapolation 
factor of 2/11 = 0.18, giving an overall background estimate for this mode of < 0.018 
event at 90%CL. 

3.3.2 Accidental backgrounds 

Accidental background was studied using data with relaxed timing requirements. 
As a consequence of the trigger system, an asymmetric timing window was used. Events 
in the time sidebands (—115 ns < At < —3 ns) and (3 ns < At < 60 ns) were used to 
extrapolate the background from the out-of-time to the in-time (—1.5 ns < At < +1.5 ns) 
signal region. At was defined as the time difference between the average track time and 
the average time of the LKr clusters forming the 7r°. Given the number of Kg decays in the 
fiducial volume and the spill length, the probability that two events overlap in time was 
~ 1 X 10^^. The out-of-time events are shown in Fig. 3. The readout system introduced a 
non-uniformity in the event time distribution which changed the effective width of the out- 
of-time window from (175 — 6) to (125 — 6) ns. Due to additional triggering effects, events in 
this window were recorded only for the 2/itight component of the trigger. An extrapolation 
factor of 1.27 was applied to take into account the relative acceptances of the 2/Xtight and 
2//ioose triggers (see Sections 3.4 and 3.5). Six events were observed in the out-of-time signal 
region, five of them in the later data-taking period when the 2/i,ioose trigger was active, 
leading to an overall accidental background estimate of ((5 x 1.27)-|-1) x 3/ (125 — 6) = 0.18 
event. 

Studies have shown that the majority of the accidental background was due to 
Ks Ti^Ti~ or Kl n^fj,^ z/^ decays which were in time with two photons from a 
K tt^tt" decay. The background for m^^ < 0.3 GeV was dominated by Kl tt^/x^ 
decays, while for 0.30 GeV < m^^u < 0.36 GeV, the background was dominated by 
Ks 7r+7r~ decays. Background where the two charged tracks came from different de- 
cays was examined by looking for events where the two tracks had the same charge. No 
same-sign event was found in the signal region, consistent with the background estimate 
above. The accidental background from associated production was found to be negligible. 

3.3.3 Background summary 

The total background was estimated to be 0.221°;!^ event in the signal region. The 
significant background contributions are summarised in Table 2. 

3.4 Trigger efficiency 

Hardware trigger: 

The trigger efficiency was dominated by the geometrical acceptance of the two-muon 
signal and by the charged vertex requirement. The efficiencies of the other trigger com- 
ponents were measured, but no correction was necessary since they were also present in 
the normalisation trigger. 

Kl —>■ 7r"'"7r~7r° events were used to measure the efficiency of the hardware coordinate 
builders and the microprocessors that reconstructed the tracks and determined the vertex 
of the decay. The efficiency of this trigger component was found to be 0.910 ± 0.001. 

The efficiency of the two-muon signal component of the trigger was estimated from 
Kl Tr+Tr"?!^ decays where both charged pions decayed to muons. The kinematics of 
these events is very similar to that of Ks — > tt^/x^/x". The efficiencies of the 2/i,tight 
and 2/xioose signals are shown in Fig. 4 as a function of m^^. A large track separation 



was effectively imposed by tlie 2/i,tight trigger requirement, resulting in a lower trigger 
efficiency at low m^u^u. As shown in Fig. 4, the dependence of the trigger efficiency on m^^ 
was consistent with that expected from Monte Carlo simulation. 

The 2yUtight and 2/iioose trigger efficiencies for Ks 7r^fi'^fi~ decays were found to 
be 0.73 ± 0.03 and 0.93 ± 0.02, respectively. Taking into account the relative amount 
of data recorded with the 2/xtight and 2fiioose triggers, and including the vertex-finder 
trigger efficiency, the overall hardware trigger efficiency correction was determined to be 
0.77 ±0.02. 

Software trigger: 

A small fraction of unbiased events was flagged by the software trigger, but not rejected. 
Using these events, the efficiency for the software trigger described in Section 2.2.4 was 
determined to be 1.00+o:o2. 

3.5 Signal Acceptance 

The signal acceptance was estimated using Monte Carlo simulation. The amplitude 
for the decay was taken to be a vector matrix element of the form [6]: 



where p, p^, p- and p+ are the four- momenta of the kaon, pion, muon and anti-muon, 
respectively, z — rnf^^/m%;, and W{z) is a form factor. As explained in [6], the W{z) 
dependence on z vanishes to lowest order and for this analysis was represented by the 
first order polynomial W{z) = as + bsz. 

The response of the LKr and HCAL detectors to muons is not well simulated by the 
Monte Carlo and the effect of the LKr and HCAL energy cuts on the muon identification 
efficiency was therefore studied using data. Prom a special run with a muon beam, the 
mean 2// identification efficiency was determined to be 0.96 ± 0.02. The dependence of 
the trigger efficiency on m^^ was also determined from data, as described in the previous 
section. 

The acceptance as a function of m^^ is shown in Fig. 5(a), before and after apply- 
ing the muon identification and trigger efficiency corrections (see Table 4). The reduced 
acceptance at low m^^ after applying these corrections is due to the 2/xtight component of 
the trigger. 

The expected m^^ distribution, taking into account the muon identification and 
trigger efficiency corrections, is shown in Fig. 5(b) for several values of bs/as, including 
a unit form factor, bs/as — 0, and bs/as — 0.4 as suggested by the Vector Meson 
Dominance (VMD) model. The reconstructed dilcpton mass distribution of the six signal 
events is also shown in Fig. 5(b) and is consistent with expectation. 

The dependence of the overall acceptance on the value of bs/as is shown in Fig. 6. 
The variation in acceptance seen in the region —5 < bs/as < —2 arises because the 
form factor develops a minimum within the kinematically allowed range of m^^^ (see 
Fig. 5(b)). The maximum variation in the acceptance over a wide range of bs/as was 
used to determine the systematic imcertainty. 

Taking into account the dependence of the trigger efficiency on m^^ and the muon 
identification efficiency, the overall acceptance was found to be 



A[Ks ttVV ] oc W{z){p + p^)''ui{p_)j^vi{p+), 



(2) 



0.081 ± 0.002{stat) ± 0.004{syst) 



(3) 



for kaons in the energy range 60 to 200 GeV decaying in the first three Ks lifetimes after 



the collimator, where a unit form factor has been assumed. For comparison, the acceptance 
before applying the trigger and muon identification corrections was 0.109 . 

3.6 Normalisation. 

The Ks flux was estimated using 132 million Kg — > tt^ti~ events from a minimum 
bias trigger. Only events that had one vertex, two tracks and no additional clusters outside 
a modified in-time window were used. Events where the charged pions had decayed were 
not rejected, and an E/p cut of less than 0.95 was used. 

The flux was corrected for tt+tt" acceptance (0.4579 ± 0.0004), software trigger effi- 
ciency (0.9989±0.0004), and the E/p cut (0.9932±0.0001), where the errors are statistical 
only. No correction was needed for the hardware trigger efficiency, for inefficiencies in the 
chambers and track reconstruction, or for losses due to multi-vertex and accidental activ- 
ity since these also applied to the signal channel. To estimate the systematic uncertainty, 
the cuts used to reject accidental activity and extra clusters were varied. A systematic 
uncertainty was assigned to take into account the 3% variation observed. 

After taking into account the prescaling factor of the trigger, the total Ks flux was 
found to be (2.50 ±0.08) x 10^° in the same fiducial volume as used for the signal channel. 

4 Result 

As shown in Fig. 1, six events were found in the signal region, with a background 
estimate of 0.22^q;^^ event. This is the first observation of the Ks — > n'^n^jji^ decay. 

The kinematic properties of the six Ks — > 7r°// candidates were consistent with 
those expected based on Monte Carlo simulation of the signal. It was also checked that the 
number of events observed when the principal analysis cuts were relaxed was consistent 
with the expected background and that there was no accumulation of events close to any 
of these cuts. 

Using the information summarised in Table 3, the Ks — > 7r°/x"'"/i~ branching ratio 
was measured to be: 

B{Ks 7rV^/i") = [2.9+iJ(stoO ± Q.2{syst)] x 10"^ (4) 

where the statistical error was obtained using the method in Ref. [7] and the systematic 
uncertainty by combining the individual errors in quadrature. 

This result is consistent within error with recent predictions based on Chiral Per- 
turbation Theory [6, 8] . 

5 Discussion 

5.1 Test of Chiral Perturbation Theory 

Chiral Perturbation Theory (ChPT) can be used to predict the branching ratio for 
Ks — > TT^t^l" and the corresponding dilepton mass spectrum, mu- The measurement 
presented here tests these predictions and constrains the parameters of the model. 

The Ks — > Ti^t^l~ branching ratios can be expressed as a function of two parame- 
ters, as and bs [6]: 

B{Ks 7r°e+e-) = [O.Ol - 0.76a5 - O.2I65 + 46.5a| + 12.90565 + 1.446|] x 10"^° (5) 
B{Ks ttV^A*") = [0.07 - 4.5205 - I.5O65 + 98.7o| + 57.70565 + 8.956|] x 10"^^ (6) 

where the total form factor is Ws{z) = GFm\{as + hsz) + Ws'^{z), z = mj^/mj^, mx is 
the kaon mass, mu is the invariant mass of the two leptons, and Ws'^{z) is expected to 



be small. Assuming VMD, which predicts bs = 0.4 05 [6], the value of the parameter \as\ 
can be obtained from the measurement of the individual Ks — > 7r°£+£~ branching ratios 
via the relations [9] 

B{Ks ^TT^e+e-) ~ 5.2xl0-^a|, (7) 
B{Ks^7r^l^+l^-) ~ 1.2xl0-^a|. (8) 

The VMD model gives an estimate of the ratio of the muon to electron branching ratios 
of 0.23 [6], while the ratio predicted with a unit form factor is 0.21 [10]. 

Within the VMD model, the value of the parameter \as\ can be obtained from the 
measurement of B{Ks — >■ 7t^ij,~^ij,~) using Eq. 8: 

This value is in agreement with the same quantity extracted from the study of Ks 
n°e+e- in Ref. [1]: 

|a5|ee = 1.061^:^? ±0.07. 

The ratio B{Ks 7i^^+^i-)/B{Ks 7r°e+e-) is found to be 0.49^2;^^ ± 0.08, in reason- 
able agreement with the VMD model prediction of 0.23. 

The 0,5 and 65 parameters can be obtained from a combined analysis of B[Ks 
TT^e^e^) and B{Kg n^fi^fi~) using Eqs. 5 and 6. The 68% confidence level contours 
in the (05,65) plane derived from the measured decay rates are shown in Fig. 7(a). The 
VMD prediction bs/as — 0.4 falls within both sets of contours. 

The constraints on as and 65 from the two branching ratio measurements have been 
combined using a maximum likelihood method, as shown in Fig. 7(b). Two regions in the 
{as, bs) plane are preferred, namely 

as = -1.6l?i, ^5= 10.81?:^ (9) 
as = 1.9tli, bs = -n.Stli. (10) 

The limited statistics do not allow a significant determination of 65 or an assessment of 
the linear dependence of the form factor on z. Table 4 shows the acceptance as a function 
of dilepton mass to facilitate the comparison of the two decay channels. 

5.2 CPV component of Kl — > 7r°/Lt+/Lt~ 

The measured branching ratio for the decay Ks 7r°/i allows the CPV con- 
tribution to the branching ratio of the corresponding K^ decay, K^ 7r°/x"'"/x~, to be 
predicted as a function of Im(Aj) to within a sign ambiguity [11]: 

B{K, ^ 7rVV-)cPV X 10^^ = Cmix ± Qnt (^^) + C^^m (^)' > (H) 
where 

Cmix = 3.1 X 10^ B{Ks ^ ttW"), CmT = 4.6 x 10^ y/B{Ks ^ n^+pi-), Cdir = 1.0. 

CiNT is the coefficient for the term due to the interference between the direct (Cdir) and 
indirect (Cmix) CPV components, and At = V'td^s- 



The predicted dependence of B{Kl 7r°/x"'"/x )cpv on Im(At) is shown in Fig. 8. 
Taking the central value of the measured branching ratio B{Ks 7r°/x+/x~) and Im(At) = 
1.36 X 10-4 [12] gives: 

B{Kl — >■ Tl^jjl^ H )cPV X 10^^ ~ 8.8niixing ± 3.3interference + l-8direct- (12) 

The predicted dependence for the T:^e^e~ channel is also shown in Fig. 8. 

The CP conserving (CPC) component of Kl Ti^i^i' decays can be constrained 
using measurements of the decay Kl — > 71^77 [13, 14]. A recent analysis based on ChPT 
obtained the prediction (5.2 ± 1.6) x 10"^^ [11]. 

Combining the CPV and the CPC components, the central value for the total 
Kl n^fi^fjr branching ratio is estimated to be 19x10"^^ or 13x 10"^^, depending on the 
sign of the interference term between the direct and the indirect (mixing) amplitudes. This 
estimate is consistent with the present experimental iippcr limit on B{Kl n^fi^fi^) of 
3.8 X 10"^° (90% CL) [15]. The corresponding central value for the total Kl 7r°e+e- 
branching ratio was estimated to be 32 x 10"^^ or 13 x 10"^^ [1]. 
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Event 


Ks momentum (GeV) 


t/ts 


rriie (GeV) 


Ks ^ ny+^i- 


1 


93.1 


0.48 


0.241 


2 


117.9 


0.71 


0.220 


3 


147.8 


1.74 


0.226 


4 


95.3 


2.49 


0.241 


5 


112.7 


0.64 


0.262 


6 


101.0 


0.96 


0.226 


Ks — > TT^e+e 


1 


84.6 


0.74 


0.291 


2 


128.2 


0.50 


0.267 


3 


114.1 


1.02 


0.173 


4 


83.9 


2.09 


0.272 


5 


130.8 


1.46 


0.303 


6 


121.2 


1.49 


0.298 


7 


94.2 


1.64 


0.253 



Table 1: Kinematics of the six events found in the signal region for Ks — > 7r°//''"/x and 
the seven events found for Ks Tr^e+e" [1]. 



Background Source 


Events 


Kl TT^TT+TT" 


"-'-0.00 


Kl /i+/i"77 




Accidentals 


o.i8;[j:l« 


Total background 


o.22ii;:i« 



Table 2: Summary of the background estimate in the signal region. 



Total Ks Flux 


(2.50 ±0.08) X 10^" 


Acceptance 


0.081 ±0.002 ±0.004 


Background 




Events 


6 



Table 3: Summary of information used to extract the Ks — > 'k'^ii'^ij, branching ratio. 



m 



(GeV) 



0.21-0.23 
0.23-0.25 
0.25-0.27 
0.27-0.29 
0.29-0.31 
0.31-0.33 
0.33-0.35 
0.35-0.37 



Acceptance 



U.U ( J_o.008 

n n79+o-oo5 

U.U ( ^_o,o06 

0.088+°;°°^ 

-0.003 
-0.003 
h0.004 
-0.006 

n nq7+o.oo3 

n 070+0.003 
U.U / J_o.012 

ni 7+0-002 

-"^ ' -0.002 



0.083 
0.084: 



K 



s 



nip 



(GeV) 



0.000-0.0365 
0.0365-0.073 
0.073-0.1095 
0.1095-0.146 
0.146-0.1825 
0.1825-0.219 
0.219-0.2555 
0.2555-0.292 
0.292-0.3285 
0.3285-0.365 



Acceptance 



0.0809 ±0.0015 
0.0877 ±0.0009 
0.0910 ±0.0007 
0.0872 ± 0.0007 
0.0785 ± 0.0006 
0.0755 ± 0.0006 
0.0703 ± 0.0007 
0.0664 ±0.0008 
0.0560 ±0.0010 
0.0430 ± 0.0020 



Table 4: Overall acceptance as a function of the dilepton mass for Ks 7r°/i+/i and 
Kg TrOg+e- [1]. 
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Figure 1: Scatter plot for the events passing all the cuts described in the text: (a) for the 
m^^ versus m^^,r plane and (b) for the m^^ versus m^^ plane. The 2.5(7 and the 6cr signal 
and control regions and the m^^ kinematic limits are also shown. 
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Figure 2: Distributions of m^^ and m^^-,^, with the cuts on m^^ and m^^^ removed, for 
data (points with error bars) and Kl 7r"'"7r~7r° Monte Carlo (histograms). The lower 
plots do not include the Ks proper hfetime cut calculated from zk- The vertical hnes 
indicate the 2.5 cr signal region. 
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Figure 3: Scatter plot of m^^ versus m^^,r for events in the out-of-time window. The 2.5a 
and the 6a signal and control regions are also shown. 




Figure 4: Trigger efficiency measured from data (points with error bars) and estimated 
from Monte Carlo (histograms) for the two-muon signal component of the trigg Gr, cis db 
function of m^^: (a) for 2/itight, and (b) for 2/iioose- 
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Figure 5: (a) Acceptance as a function of m^^, including (hatched histogram) and not 
including (open histogram) the trigger efficiency; (b) Expected m^^ distributions of ac- 
cepted events for several values of bs/as after taking into account the trigger efficiency. 
The m^^ distribution of the six signal events is also shown. 
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Figure 6: Overall acceptance as a function of bs/as for the Ks vr'^yU+yU" (upper curve) 
and Ks —>■ vr^e^e" (lower curve) channels. 
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Figure 7: (a) Allowed regions of and 65 determined from the observed number 
of Ks 7r*'/i+yU~ and Ks —>■ Ti^e^e" events separately. The region between the inner 
and outer solid (dashed) elliptical contours is the allowed region for Ks vr^/i+yU" 
{Ks — > 7i^e~^e~) at 68% CL. (b) Allowed regions of as and bs for the Ks 7r°/i+/i~ 
and Ks — > TT^e^e~ channels combined. The inner (outer) contour of each pair delimits 
the Icr (2(t) allowed region from the combined log-likelihood. The dashed straight line in 
both plots corresponds to 65 = 0.405, as predicted by the VMD model. 
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Figure 8: Predicted CPV component of the K^ vr^/i+yU" (solid curves) and 
Kl 7r°e"^e" (dashed curves) branching ratios as a function of Im(Aj) assuming (a) 
as < and (b) as > 0. Each pair of curves delimits the allowed range derived from the 
±1(7 measured values of \as\- The vertical shaded band shows the world average value of 
Im(A() with its uncertainty [12]. 



